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In this article, we introduce multiple turbo codes and a suitable decoder structure 
derived from an approximation to the maximum a posteriori probability (MAP) 
decision rule, which is substantially different from the decoder for two-code- based 
encoders. We analyze the effect of interleaver choice on the weight distribution of 
the code, and we describe simulation results on the improved performance of these 
new codes. 


I. Introduction 

Coding theorists have traditionally attacked the problem of designing good codes by developing codes 
with a lot of structure, which lends itself to feasible decoders, although coding theory suggests that codes 
chosen “at random” should perform well if their block size is large enough. The challenge to find practical 
decoders for “almost” random, large codes has not been seriously considered until recently. Perhaps the 
most exciting and potentially important development in coding theory in recent years has been the 
dramatic announcement of “turbo codes” by Berrou et al. in 1993 [1], The announced performance of 
these codes was so good that the initial reaction of the coding establishment was deep skepticism, but 
recently researchers around the world have been able to reproduce those results [3,4], The introduction 
of turbo codes has opened a whole new way of looking at the problem of constructing good codes and 
decoding them with low complexity. 

It is claimed these codes achieve near-Shannon-limit error correction performance with relatively simple 
component codes and large interleavers. A required E b /N 0 of 0.7 dB was reported for a bit error rate 
(BER) of 10 [1], However, some important details that are necessary to reproduce these results were 

omitted. The purpose of this article is to shed some light on the accuracy of these claims and to extend 
these results to multiple turbo codes with more than two component codes. 

The original turbo decoder scheme, for two component codes, operates in serial mode. For multiple- 
code turbo codes, we found that the decoder, based on the optimum maximum a posteriori (MAP) rule, 
must operate in parallel mode, and we derived the appropriate metric, as illustrated in Section III. 


II. Parallel Concatenation of Convolutional Codes 

The codes considered in this article consist of the parallel concatenation of multiple convolutional 
codes with random interleavers (permutations) at the input of each encoder. This extends the analysis 
reported in [4], which considered turbo codes formed from just two constituent codes. Figure 1 illustrates 
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Fig. 1. Example of encoder with three codes. 


a particular example that will be used in this article to verify the performance of these codes. The 
encoder contains three recursive binary convolutional encoders, with Mi, M 2 and M 3 memory cells, 
respectively. In general, the three component encoders may not be identical and may not have identical 
code rates. The first component encoder operates directly (or through 7Ti) on the information bit sequence 
u — ( w 1j * * * i u n) of length N , producing the two output sequences x iz and x lp . The second component 
encoder operates on a reordered sequence of information bits, u 2 , produced by an interleaver, 7 r 2 , of 
length TV , and outputs the sequence x 2p . Similarly, subsequent component encoders operate on a reordered 
sequence of information bits, u,, produced by interleaver tt J 5 and output the sequence x jp . The interleaver 
is a pseudorandom block scrambler defined by a permutation of N elements with no repetitions: A 
complete block is read into the the interleaver and read out in a specified (fixed) random order. The 
same interleaver is used repeatedly for all subsequent blocks. Figure 1 shows an example where a rate 
r = \/n = 1/4 code is generated by three component codes with Mi = M 2 = M 3 = M = 2, producing 
the outputs xij = u, xi p = u • gb/g a , x 2p = u 2 • gb/g a , and x 3p = u 3 • gb/g a (here tt\ is assumed to be 
an identity, i.e., no permutation), where the generator polynomials g a and g^ have octal representation 
{7) octal and ( 5) octal, respectively. Note that various code rates can be obtained by proper puncturing 
°f x ipT x 2 p ? x 3 pi and even xj t if the decoder works (for an example, see Section IV). The design of the 
constituent convolutional codes, which are not necessarily optimum convolutional codes, is still under 
investigation. It was suggested in [5] that good codes are obtained if g a is a primitive polynomial. 

We use the encoder in Fig. 1 to generate an (n(N + M), TV) block code, where the M tail bits of 
code 2 and code 3 are not transmitted. Since the component encoders are recursive, it is not sufficient to 
set the last M information bits to zero in order to drive the encoder to the all-zero state, i.e., to terminate 
the trellis. The termination (tail) sequence depends on the state of each component encoder after N bits, 
which makes it impossible to terminate all component encoders with M predetermined tail bits. This 
issue, which had not been resolved in previously proposed turbo code implementations, can be dealt with 
by applying the method described in [4], which is valid for any number of component codes. 
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A. Weight Distribution 

In order to estimate the performance of a code, it is necessary to have information about its minimum 
distance, weight distribution, or actual code geometry, depending on the accuracy required for the bounds 
or approximations. The challenge is in finding the pairing of codewords from each individual encoder, 
induced by a particular set of interleavers. Intuitively, we would like to avoid joining low-weight codewords 
from one encoder with low- weight words from the other encoders. In the example of Fig. 1, the component 
codes have minimum distances 5, 2, and 2. This will produce a worst-case minimum distance of 9 for the 
overall code. Note that this would be unavoidable if the encoders were not recursive since, in this case, the 
minimum weight word for all three encoders is generated by the input sequence u = (00 ■ 0000100 ■ ■ • 000) 

with a single “1,” which will appear again in the other encoders, for any choice of interleavers. This 
motivates the use of recursive encoders, where the key ingredient is the recursiveness and not the fact 
that the encoders are systematic. For our example, the input sequence u = (00 ••• 00100100 • • • 000) 
generates a low-weight codeword with weight 6 for the first encoder. If the interleavers do not “break ’ 
this input pattern, the resulting codeword’s weight will be 14. In general, weight-2 sequences with 
2 4- 3t zeros separating the l’s would result in a total weight of 14 + 6 1 if there were no permutations. By 
contrast, if the number of zeros between the ones is not of this form, the encoded output is nonterminating 
until the end of the block, and its encoded weight is very large unless the sequence occurs near the end 
of the block. 

With permutations before the second and third encoders, a weight-2 sequence with its l’s separated 
by 2 + 3f i zeros will be permuted into two other weight-2 sequences with l’s separated by 2 + 3t, zeros, 
i = 2, 3, where each t, is defined as a multiple of 1/3. If any U is not an integer, the corresponding encoded 
output will have a high weight because then the convolutional code output is nonterminating (until the 
end of the block). If all U's are integers, the total encoded weight will be 14 + 2 £ i=1 U. Thus, one of the 
considerations in designing the interleaver is to avoid integer triplets (ti T *2 > ^3) are simultaneously 

small in all three components. In fact, it would be nice to design an interleaver to guarantee that the 
smallest value of £< =1 ( for inte S er U) grows with the block size N. 

For comparison, we consider the same encoder structure in Fig. 1, except with the roles of g a and 
g b reversed. Now the minimum distances of the three component codes are 5, 3, and 3, producing an 
overall minimum distance of 11 for the total code without any permutations. This is apparently a better 
code, but it turns out to be inferior as a turbo code. This paradox is explained by again considering 
the critical weight-2 data sequences. For this code, weight-2 sequences with 1 + 2ti zeros separating the 
two l’s produce self-terminating output and, hence, low- weight encoded words. In the turbo encoder, 
such sequences will be permuted to have separations 1 + 2t„i = 2,3, for the second and third encoders, 
where now each t x is defined as a multiple of 1/2. But now the total encoded weight for integer triplets 
(tut 2 ,t 3 ) is 11 + t i- N o tice how this wei S ht g rows onl V half 38 fa8t with £*=1 38 the Previously 
calculated weight for the original code. If £- =1 U can be made to grow with block size by the proper 
choice of an interleaver, then clearly it is important to choose component codes that cause the overall 
weight to grow as fast as possible with the individual separations t x . This consideration outweighs the 
criterion of selecting component codes that would produce the highest minimum distance if unpermuted. 

There are also many weight-n, n = 3, 4, 5, • - -, data sequences that produce self-terminating output 
and, hence, low encoded weight. However, as argued below, these sequences are much more likely to be 
broken up by the random interleavers than the weight-2 sequences and are, therefore, likely to produce 
nonterminating output from at least one of the encoders. Thus, turbo code structures that would have 
low minimum distances if unpermuted can still perform well if the low-weight codewords of the component 
codes are produced by input sequences with weight higher than two. 

B. Random Interleavers 

Now we briefly examine the issue of whether one or more random interleavers can avoid matching small 
separations between the l’s of a weight-2 data sequence with equally small separations between the 1 s of 
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ltspermuted version (s)- Consider, for example, a particular weight-2 data sequence (• • ■ 001001000 • • ■) 
which corresponds to a low- weight codeword in each of the encoders of Fig. 1. If we randomly select an 
mteHeave, of size Al, .he probability that .his sequence will be pernru.e/iuto another TJ/i "f ,h. 
e form is roughly 2/N (assuming that TV is large and ignoring minor edge effects). The probability 

* leaSt «» o.iginal sequence 


•001001000- • •) within the block size of TV is approximately 1 - (1 - 2/N) 


N 


,i , i . . . . - - rr . V aw.. i ^v V x — &/ iy j ~1 — e mis implies 

that the minimum distance of a two-code turbo code constructed with a random permutation is not likely 
to be much higher than the encoded weight of such an unpermuted weight-2 data sequence, e.g., 14 for the 
code m Fig. 1 (For the worst-case permutations, the d mtn of the code is still 9, but these permutations 
, G hlg !\ Unikey lf chosen randomly.) By contrast, if we use three codes and two different interleavers 
17K a" Particular sequence (• • • 001001000 • • •) will be reproduced by both interleavers is 
y (2/N) . Now the probability of finding such an unfortunate data sequence somewhere within the 
block of size TV is roughly 1 - [l - (2/TV) 2 ]" » 4/TV. Thus, it is probable that a three-code turbo code 
using two random interleavers will see an increase in its minimum distance beyond the encoded weight 
o an unpermuted weight-2 data sequence. This argument can be extended to account for other weight-2 
data sequences that may also produce low-weight codewords, e.g., (• • ■ 00100(000)* 1000 • • •), for the code 
in lg . 1 . For comparison, let us consider a weight-3 data sequence such as (• • • 0011100 • ■ ■), which for our 
example corresponds to the minimum distance of the code (using no permutations). The probability that 
is sequence is reproduced with one random interleaver is roughly 6//V 2 , and the probability that some 
sequence of the form (• • • 0011100 • • •) is paired with another of the same form is 1 - (1 -6/N 2 ) n k 6/N. 
lius, for large block sizes, the bad weight-3 data sequences have a small probability of being matched with 
ad weight- 3 permuted data sequences, even in a two-code system. For a turbo code using three codes and 
two random interleavers, this probability is even smaller, 1 - [l - (6/N 2 ) 2 ] N « 36/N 3 This imnlies that 
the minimum distance codeword of the turbo code in Fig. 1 i mom likely to result from a w“g£j data 
sequence of the form (• • ■ 001001000 • • •) than from the weight-3 sequence (■ ■ ■ 001 1 100 • • •) that produces 
the minimum distance in the unpermuted version of the same code. Higher weight sequences have an 
even smaller probability of reproducing themselves after being passed through the random interleavers. 

For a turbo code using q codes and q- 1 interleavers, the probability that a weight-n data sequence will 
be reproduced somewhere within the block by all g -1 permutations is of the form 1-Tl - (P/N n ~ l )i~ l ] N 
where (3 is a number that depends on the weight-n data sequence but does not increase with block size 
IV. For large TV, this probability is proportional to (1/2 V)"«-»-*, which falls off rapidly with TV when n 
and q are greater than two. Furthermore, the symmetry of this expression indicates that increasing either 

probIbnitv° f th<3 data S6qUenCe ” ° r tHe number of codes V has roughly the same effect on lowering this 

In summary, from the above arguments, we conclude that weight-2 data sequences are an important 
factor in the design of the component codes, and that higher weight sequences have successively decreasing 
importance. Also increasing the number of codes and, correspondingly, the number of interleavers, makes 
i more and more likely that the bad input sequences will be broken up by one or more of the permutations. 

The minimum distance is not the most important characteristic of the turbo code, except for its 
asymptotic performance, at very high E b /N 0 . At moderate signal- to- noise ratios (SNRs), the weight 
distribution for the first several possible weights is necessary to compute the code performance. Estimating 
e complete weight distribution of these codes for large TV and fixed interleavers is still an open problem 
However, it is possible to estimate the weight distribution for large TV for random interleavers by using 
probabilistic arguments. (See [4] for further considerations on the weight distribution). 

C. Design of Nonrandom and Partially Random Interleavers 

Interleavers should be capable of spreading low-weight input sequences so that the resulting codeword 
as high weight. Block interleavers, defined by a matrix with v r rows and u c columns such that N = u r xu 
may fail to spread certain sequences. For example, the weight-4 sequence shown in Fig. 2 cannot be broken 


-2 


This implies 
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Fig. 2. Example where a block Interleaver fails to 
“break” the input sequence. 

bv a block interleaver. In order to break such sequences, random interleavers a.re desirable, as discussed 
above (A method for the design of nonrandom interleavers is discussed in [3]) Block interleavers are 
effective if the low-weight sequence is confined to a row. If low-weight sequences (which can be regarded as 
the combination of lower-weight sequences) are confined to several consecutive rows then the columns 
of the interleaver should be sent in a specified order to spread as much as possible the low-weight sequence^ 
A method for reordering the columns is given in [7]. This method guarantees that for any number of 
columns v c = a<j + r, (r < a - 1), the minimum separation between data entries is q - 1, where a is 
the number of columns affected by a burst. However, as can be observed in the example in lg. , e 

sequence 1001 will still appear at the input of the encoders for any possible column permutation. Only 
if we permute the rows of the interleaver in addition to its columns is it possible to break the low-weight 
sequences The method in [7] can be used again for the permutation of rows Appropriate selection of 
and q for rows and columns depends on the particular set of codes used and on the specific low-weight 
sequences that we would like to break. 

We have also designed semirandom permutations (interleavers) by generating random integers t, 
, <i<N, without replacement. We define an “5-random" permutation as follows: Each random* 
selected integer is compared to S previously selected integers. If the current refection is equa Uomry 
5 previous selections within a distance of ±S , then the current selection is rejected. I his process is 
repeated until all N integers are selected. The searching time for this i algorithm increases w|h 5 and 
is not guaranteed to finish successfully. However, we have observed that choosing S < y/N/2 usually 
produces a solution in a reasonable time. Note that for 5 = 1, we have a purely random interleaver. In 
the simulations, we used S = 31 with block size N — 4096. 


III. Turbo Decoding for Multiple Codes 


In this section, we consider decoding algorithms for multiple-code turbo codes. In the ad- 

vantage of using three or more constituent codes is that the corresponding two or more interleavers have 
TSer chance to break sequences that were not broken by another interleaver. The disadvan age is 
that for an overall desired code rate, each code must be punctured more, resulting in weaker constituent 
codes. In our experiments, we have used randomly selected interleavers and interleavers based on the 

row— column permutation described above. 

A. Turbo Decoding Configurations 

The turbo decoding configuration proposed in [1] for two codes is shown schematically in Fig. 3. This 
configuration operates in sefial mode, i.e., “Dec 1” processes data ^before “Dec 

and so on An obvious extension of this configuration to three codes is shown in Fig. 4(a), which also 
operates in serial mode. But, with more than two codes, there are other possible configurations, such as 
that shown in Fig. 4(b), where “Dec 1” communicates with the other decoders, but these decoders 
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Fig. 3. Decoding structure for two codes. 


(a) 



H9 rat eoiiLT/Ki* decodin 9 Structures for three codes: 

(a) serial, (b) master and slave, and (c) parallel. 

not exchange information between each other This “mast-nr .,,,,1 e i„ » 

serial-parallel mode, since all other decndpr* , ,, c d configuration operates in a mixed 

shown in Fig. 4(c) is that all decoders > f CXCep ® ^ rst operate in parallel. Another possibility, 

allowed in these structures since they caus^deirrada^ ° glVen tlme ‘ Note that self loo P s are not 

feedback). We are not ? deCOding pr ° CeSS ( positive 

better? Our selection of the b^cSSn hs . Wh '' h “"juration performs 

analysis of the minimum-bit-error decoding rule (MAP algSm), as^e^'belo^** 1 * ““ 


B. Turbo Decoding for Multiple Codes 


bits 1 : = {o ; <* ,„f„,ma,io„ 

symbols y: ' alg ° n,h '” |6 ' prov ‘ des the >°S likelihood ratio L k , given the received 


■P(nic = l|y) 
P(u k = 


( 1 ) 


= 1^ gu= ut= iP(y|u)n^ fc P(uj) 
Eu:u t =0^(y|u) I] j+kPfrj) 


+ log 


P{ Uk = 1 ) 
P(u k = 0) 


( 2 ) 
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y 1( =p( 2u - 1 ) + n i; 

y 1p =p(2x 1p" 1) + n 1p 


„ f . fVn , 0 ) when the a priori probabilities P(uj) are nonuniform, the modified 
Sf" Jus, than the version considered in W . Therefore, in this a, tide, we 

use the modified MAP algorithm of [2], as we did in [4], 

The channel model Is shown in Fig. 5, where the n.e's and the »*. .re 
distributed (i.i.d.) sero-mean Gaussian random variabte wdh un« var.a “restrict 
SNR. The same model is used for each encode . Tb expton the to d g P the 

ourselves to three codes, but extension to severai cf.des ^s rmght oT outputs 

notation, consider the combination of permuter optimum bit 

V ,-012 3(xo = U) and the corresponding received sequences y t , t ’ 

to, siM metri'c on each Lit is (for data with uniform a prior, probabilities) 


y , P(vnlu)P(yi|u) J> (y 2l u ) p (y3l u ) 
Lk = og Euu k =o- p (yoi u ) p ( yi|u)P(y2,u)p(y3 ' u) 


( 3 ) 


a. p c\\ f or large N because the permutations w 2 , ^3 imply that y 2 
but in practice, we cannot comput q. ( ) SuDDose that we evaluate P(y,|u), , = 0,2, 3 

and y 3 are no longer simple convolutional encodings of u. Suppose 

in Eq. (3) using Bayes’ rule and using the following approximation. 


N 


P( u|yt) » n Pi(y-k) 


( 4 ) 


k= 1 


Note that P(u|y.) is not separable in general. 

holds with equality. If such an approximation ^ j’ A reasona ble criterion for this approximation 

i = 2 and z = 3 (by Bayes’ rule) to complete the algonthim A reas Define by 

is to choose TIL, A(«fc) such that it minimizes the Kullback distance 


Pi(u k ) = 


e UkL t k 

1 4 - e^ ,k 


where u k € { 0 , !}• Then the Kullback distance is given by 


( 5 ) 


F(L*) = E 


u k L, k 

p 4 —^ = 1 


nti(l + e^) l0g nr=i(! + e^)P(n\y t ) 


V N ti k L. 

p/-Jk = l 


( 6 ) 
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Minimizing F(L,) involves forward and backward recursions analogous to the MAP decoding algorithm 
but we have not attempted this approach in this work. Instead of using Eq. (6) to obtL {P } or’ 
equivalently, we use Eqs. (4) and (5) for i = 0,2,3 (by Bayes’ rule) to express Eq. ( 3 )“ 1 ^ ’ 


Lk /(yi i Lq, i< 2 , L 3 , k) + Lg k + L/ 2 k 4- Li 3 k 


(7) 


where Lo k = 2 py 0 fc and 


/(yi.I'O, L 2 , L 3 , k) = log 


Eu:v t =i p (yi|u) rijvfc e u i( L °i+ L 2i+ L 3i) 
52u-.u k =o p (yi|u) EU e“Aioj+Ijj+I 3 j) 


We can use Eqs. (4) and (5) again, but this time for i = 0, 1,3, to express Eq. (3) as 


(8) 


T'fc — /(y2> Lq, Li, £3, fc) + Z ofc + Li* + L 3k 


(9) 


and similarly, 


/(ys.LcLj^a,^) + Z 0fc + Z lfc + Z 2fc 


( 10 ) 


A solution to Eqs. (7), (9), and (10) is 

£ifc = /(yi,Lo,L 2l L3,*); =/(y a ,Lo l L lj L 3 ,*); L 3 * = /(y 3 ,Lo,Li,L 2 ,ifc) ( 11 ) 

for k : = 1,2,- -.AT, provided that a solution to Eq. (11) does indeed exist. The final decision is then 


L k = Lok + L ik 4 - Z 2 fc 4- L 3 k 


( 12 ) 


which is passed through a hard limiter with zero threshold. We attempted 
in Eq. (11) for Li, L 2 , and L 3 by using the iterative procedure 


to soive me nonlinear equations 


f(m + l) _ (m) rt f f (m) f (m) ,, 

^1 k - a i j{y 1,^0, \L k 3 ; ,k) 

for k = 1 2, - , AT, iterating on m. Similar recursions hold for L { ™ ] and L { ™\ The gain a\ m) should 
be equal to one but we noticed experimentally that better convergence can be obtained by optimizing 
this gam for each iteration, starting from a value slightly less than one and increasing toward one with 
the iterations, as is often done in simulated annealing methods. We start the recursion with the initial 


(13) 


nditmn Lj - L 2 - L 3 - L 0 . For the computation of /(•), we use the modified MAP algorithm 
described in [4] with permuters (direct and inverse) where needed, as shown in Fig. 6 for block 
decoder 2. The MAP algorithm always starts and ends at the all-zero state since we always terminate 
the trellis as described in [4] Similar structures apply for block decoder 1 (we assumed 7 rj= I identity 
however, any tt, can be used) and block decoder 3. The overall decoder is composed of block decoders 


1 Note that the components of the L*’s corresponding to the tail bits, i.e., L ik , for k = N + 1 
for all iterations. ’ 


• 1 N 4 A/, are set to zero 
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Fig. 6. Structure of block decoder 2. 


. , in Fig 4(c) wh ich can be implemented as a pipeline or by feedback. We proposed an 

connected as in g. I )> n Mnl At this point further approximation for turbo decoding 

Si' ti, replacing with f °” = »• 1 A 

can be used to, L» and L zk in Eq. (11). This suboptimum d^er then spends to a turbo deco 
that uses soft output Viterbi (SOVA)-type decoders rather than MAP decoders. 


C. Multiple-Code Algorithm Applied to Two Codes 

For turbo codes with only two constituent codes, Eq. (13) reduces to 


f (m+1) 
L \ k 


a 


(m) 


/(yi,L 0 ,L^ m) ,fe) 


4 ? +1) = * 


(T " ) /(y 2 ,L 0 ,L ( r ) ,fc) 


r r. 19 N and m — 1 2 • • • where, for each iteration, a< m) and a< m) can be optimized (simulated 
structure reduces to duplicate copies of the structure in Fig. 3 (i.e. , to the m ) • 



Fig. 7. Parallel structure for two codes. 

If we optimize and our method for 

in (1), which requires estimates of the “nances of * was subtracted from In. 

errors. In the method proposed in [2] the eceiv^ syste^ ^ [2] ^ ^ but the rec eived 

which results in performance degra • 11 decoder 2 In [41 we argued that there is no 

“systematic” observation was interleaved and provided to decoder 2. in [4|, we arg 
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need to interleave the received “systematic” observation and provide it to decoder 2 since L», does this 
job. It seems that our proposed method with nnr i (m) , 5 °* c * oes ^ 1S 

same performance repoLd Tp]"m ,l l/2 “dee ’ eqUi *' ‘° ' " ““ “ d * d ‘ i ™ ** 

D. Terminated Parallel Convolutional Codes as Block Codes 

Consider the combination of permuter and encoder as a linear block rode Define p * f u 
mat™ of the terminated convolutional code i. Then the overall generator matrix for thr« parallel coto 


G=[I 71-jP, TT 2 P 2 TTgPg] 


^Z‘r l th 7Tr“ i0nS ( “ ter,eav ^). In order to maximize the minimum distance of the code 


IV. Performance and Simulation Results 


For comparison with the new results on three-code turbo codes, we reproduce in Fig 8 the performance 
obtained m d b, using twmcode K = 5 turbo codes with generators (l.Wsa) 

“ d and with random permutations of lengths N = 4096 and iv = 16384 The brat 

ren!!t7 nCe “ F .‘ g 'i ™ a PP roximat ely 0.7 dB from the Shannon limit at BER = l(r 4 We also 

Zt STJS r e r" K , °‘“T' W ^ — « -»l- wi* b „„ec,ual raiwlth 
9 = (25) To show th ( 1 , 9t/9s.9 e /9.) and ( 94 / 9 .), where g. = (3T)„ to ,, g b = (33)„ a „,, and 
^ 11 r 1 /o u th ^ lt 1S possibIe not t0 send uncoded information for both codes we used an 

everali rate 1/2 turbo code using two codes with K = 2 (differential encoder) with genemter (a!/, ) 

- /-ir, ° A Ct u- a,K 9b ~ ^ octal ’ and a K ~ 5 co de with generator {g b /g a ), where q = ( 23 ) and 

gb (33 )octai- A bit error rate of 10 5 was achieved at BSNR n Sfi wr • c f octal 

of length N = 16,384 with 5 = 40. " °' 85 dB USmg an S - random permutation 

A. Three Codes 

tUrb ° COd “ With random is shown in Fig. 9 

-ru 4 ^ 96 ' 7 h fi C ° de Uses three recursi ve codes shown in Fig. 1 with constraint length K — 
The second code uses three recursive codes with K = 4 a = (131 , ! ' ^ f " 3 ' 

W A / See “| n Fi ?- 9 ’ the performance of the K = 4 code was improved by going from 20 to 30 iterations 
e ound that the performance could also be improved by using an 5 -random interleaver with 5 = 31. ’ 


V. Conclusions 

uai^T Sh ° W ” h °” , thrCe - code turb » »*» and decoders can be used to further improve the coding 
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Previous estimates on the degradations due to a subcarrier loop assume a square- 
wave subcarrier . This article provides a closed- form expression for the degradations 
due to the subcarrier loop when a finite number of harmonics are used to demod- 
ulate the subcarrier, as in the case of the buffered telemetry demodulator. We 
compared the degradations using a square wave and using finite harmonics in the 
subcarrier demodulation and found that, for a low loop signal-to-noise ratio, using 
finite harmonics leads to a lower degradation. The analysis is under the assumption 
that the phase noise in the subcarrier (SC) loop has a Tikhonov distribution. This 
assumption is valid for first-order loops. 


I. Introduction 

In an imperfect subcarrier demodulation, the difference between the phase of the reference signal and 
that of the subcarrier of the received signal causes the signal power to degrade while the noise power 
remains the same. This degradation is measured as the ratio of the reduced symbol energy-to- noise density 
ratio (E s /No), or symbol signal-to-noise ratio (SNR), to the symbol SNR of an ideal demodulation where 
the phase difference is zero. The degradations due to the subcarrier loop were previously computed 
assuming a square wave [3]. This assumption is inappropriate in the case where only a finite number 
of harmonics of the subcarrier are there to be demodulated, as in the buffered telemetry demodulator 
(BTD) [2]. This article provides a closed-form expression for computing the degradation due to a finite- 
harmonic subcarrier tracking loop. Numerically, we found that, for low loop SNR cases, we actually have 
less degradation using a finite number of harmonics than using “all” the harmonics, namely, the square 
wave. The degradation due solely to the subcarrier loop using four harmonics is 0.15 to 0.3 dB lower 
than that using a square wave for loop SNRs in the range of 14 to 30 dB. 

At first glance, the above may seem to contradict the intuition that the more harmonics we use, the 
higher the SNR we should get. This intuition is correct when the loop SNR is high, that is, when the jitter 
of the phase difference (between the true and the reference phases) is low. At low loop SNRs, however, 
we have a different scenario. 

To explain this, let us first take a look at how the subcarriers are demodulated. A square- wave 
subcarrier is demodulated by multiplying the received signal by a square-wave reference signal. When we 
only have a finite number of harmonics of the square-wave subcarrier, the current design for the BTD [2] 
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